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We study the possible formation of large (mesoscopic) molecular ions in an ultracold degenerate 
bosonic gas doped with charged partides (ions). We show that the polarization potentials produced 
by the iònic impurities are capable of capturing hundreds of àtoms into loosely bound states. We 
describe the spontaneous formation of these hollow molecular ions via phonon emission and suggest 
an optical technique for coherent stimulated transitions of free àtoms into a specific bound state. 
These results open up new interesting possibilities for manipulating tightly confined ensembles. 

PACS numbers(s): 33.90.+h, 34.10. +x, 36.40.-c, 51.90.+r 



The experimental realization of Bose-Einstein conden- 
sation (BEC) in atòmic samples ji] have stimulated a new 
wave of theoretical and experimental studies of degener- 
ate systems in the dilute and weakly interacting regimc. 
A broad range of techniques from atòmic physics and òp- 
tics have allowed the accurate manipulation of ultracold 
samples of neutral àtoms. In particular, various proper- 
ties of degenerate gases have been explored, such as their 
excitation modes Q, superfluidity phenomena J^Q, and 
controlled vortex creation || . Recent studies are probing 
ultracold atòmic systems in which elèctric charges may 
play an important role || . These include ultracold plas- 
mas [Q , ultracold Rydberg gases Q , as well as ionization 
experiments in BEC |J. 

In this Letter, we explore theoretically the behavior 
of a dilute atòmic Bose-Einstein condensate doped with 
iònic impurities. We show that the polarization inter- 
action in a condensate can lead to the capture of large 
numbers of àtoms into weakly bound states, resulting 
in the ràpid formation of mesoscopically large molecular 
ions. We study the spontaneous dynamics of the molecu- 
lar ion formation and show that the degenerate nature of 
the condensate, as well as the properties of collective ex- 
citations (phonons), play an important role. We furthcr 
describe a coherent optical technique to prepare molecu- 
lar ions in specific states. Beside the fundamental interest 
of studying the formation of such large many-body ob- 
jects, the effects described here may open up new ways 
to manipulate cold àtoms. In particular, the charged and 
tightly trapped atòmic cloud represents a microtrap that 
could easily be manipulated and "transported" by exter- 
nal fields. Controlled mechanisms to manipulate tightly 
confined, strongly interacting àtoms may also allow for 
new approaches for quantum information processing and 
for studies of other fascinating phenomena such as quan- 
tum phase transitions. Before proceeding we note an 
interesting analogy to early studies involving charged im- 
purities in superfluid helium, where electron bubbles and 
ion "snowballs" were predicted and observed |ÏÏJ . 

We considcr the situation in which fcw ions with ul- 



tralow kinetic energy are introduced into an atòmic BEC. 
This can be realized, e.g., by rapidly ionizing àtoms from 
the condensate using làsers in a process where the ejected 
electrons carry essentially all the kinetic energy. This 
would leave the BEC doped with few ions, and the BEC- 
ion system in a non-equilibrium state. Alternatively, it 
could be possible to introduce charged impurities via con- 
trolled processes involving either a combination of ion 
and atom traps, or using surface traps on semiconduc- 
tor surfaces. As a relaxation process, large numbers of 
àtoms from the condensate can be captured into loosely 
bound states of the polarization potentials, rapidly form- 
ing shells of àtoms around ions. Such a process occurs 
spontaneously through collisions of condensate àtoms in 
which àtoms are stimulated down into the molecular ion 
bound state, and the excess energy is carried away by the 
condensate collective excitations. 

We are interested in the limit T — > 0, and for simplic- 
ity, we consider a homogeneous BEC with the neutral gas 
being the parent atom of the doping ion [ pd| . The ion will 
polarize a nearby atom (separated by a distance r) and 
the two will interact via a polarization potential behav- 
ing asymptotically as —0^/2^, where C4 is the dipole 
polarizability of the neutral atom. 

Under the conditions mentioned above, two types of 
scattering processes involving the ion are possible: elàstic 
and super-elastic collisions. The energy of the colliding 
partners is not changed in the first type, while kinetic 
energy is released in the second type. Elàstic scattering 
of an atom and an ion is often described as the contribu- 
tions from two separate processes ||[Ï2| : resonant charge 
transfer scattering (the electron changes center after the 
collision), and pure elàstic scattering (the electron rc- 
mains with the same center). According to Wigner's 
threshold laws ]Ï3| ], the rate of all elàstic processes van- 
ishes as T — > 0. The second type of scattering process, 
i.e. super-elastic scattering, corresponds to an inelastic 
process where one condensate atom is captured by the 
ion and where kinetic energy is released via the emission 
of phonons. Here, condensate àtoms are accelcratcd by 
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the iònic field (see Fig. ||a), and after the collision, onc 
of the neutral àtoms is captured by the polarization po- 
tential, and the kinetic energy released is shared by the 
remaining free condensate àtoms and the newly formed 
molecular ion |y| . Inelastic collisions with excitations of 
BEC collective modes usually cannot be decribed by sim- 
ple binary collisions, unless in the limit where the binding 
energy of the uppermost bound level of the molecular ion 
is much larger than the chemical potential of the BEC. 
According to Wigner's threshold laws JÏ3|], the rate for 
super-elastic processes tends to a constant at zero tem- 
perature, and this scattering process will be dominant at 
T — ► 0. Contrary to the slow-down of impurities in con- 
densates, where phonon radiation is not possible below 
some critical velocity (because of momentum-energy con- 
servation the capture of àtoms by the long range po- 
tential created by the ion corresponds to free-bound tran- 
sitions, and does not suffer from this restriction: phonon 
emission takes place at any velocity. 

We describe the evolution of the number of àtoms N v 
in the bound level v of the polarization potential by a 
kinetic equation (see Fig. |^a) 



dN v 



W™ P (N V + 1) - (W, down + W^)N V , (1) 



where WJ ap is the capture rate from the condensate, and 
Wy° wn and W^ p are the loss rates to more deeply bound 
states and back to the condensate, respectively. The fac- 
tor N v + 1 comes from the bosonic nature of the àtoms: 
the capture will be bose-enhanced (or stimulated), as op- 
posed to the depletion of the level (proportional to N v ). 

Let us first calculate W^ ap . As mentioned above, the 
capture rate is equal to the phonon emission rate. For Bo- 
goliobov quasi-particles (or phonons), the emission rate 
of phonons with momemtum q and energj 
estimated using the Fermi Golden rule 0JÏ 
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where V is the volume occupied by the condensate, 
fj, c = Antena/m its chemical potential (n, m, and a are 
the density, mass, and scattering length of the condensate 
àtoms, respectively), and riq the phonon occupation num- 
ber. Herc, I(q) is the square of the form factor for transi- 
tions from a continuum to a bound state. If we represent 
the initial condensate state by the single particle wave- 
function í'o (^0 an d the final bound state by the wavefunc- 
tion í' t ,(r), with energy difference As — Sq — e v — huj^, 
then 



I(q) = N 



d 3 r*:(f)e- l? "*o(r) 



(3) 



where N — nV is the number of àtoms in the condensate. 

At zero-temperature in an infinite homogeneous sys- 
tem, we set e = (k = 0), and we have pi 
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where the binding energy of the final bound state is 
related to the extent of the s-wavefunction ^> v (f) via 
e v = —h 2 /2/ia 2 7 with [i being the system reduced mass. 
Note that we have Ae = h 2 /2/ia 2 . In an isotropic Bose 
gas, I{q) depends only on the magnitude of q, and with 
Eq.@ í m=8na 3 v n/(l + q 2 a 2 v ) 2 . 

The capture rate is obtained by integrating over all 
possible phonon states, i.e. J d 3 q w e mis.(<7) ^ / /(27r) 3 , and 



wr = iVn f^ )3/2 Ko + Dl fa) , (5) 

h V a* a v v 1 + í 

where Tiui qo = Ae, and £ = /i c /Ae = 8irna 2 a[j,/m. For 
Bogoliobov phonons, TiLj q = Tiqs^fl + (hq/2ms) 2 , and we 
can express qo as 



<lla 2 v ^ = ^ÏTe-i, 



(6) 



where the sound velocity s is related to the chemical po- 
tential by /x c = ms 2 , and W^ v becomes 
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For infinitly massive ions, we have m/n = 1 p7| . Note 
that the two limiting vàlues of £ correspond to two dif- 
ferent physical conditions: £ — > oo implies qo — > 0, i.e. a 
phonon-like regime, and £ — > implies dilute conditions 
with a binary collision regime. In Fig. ||, we illustrate 
W^ p as a function of £, and show its behavior for the 
asymptotic regimes: WJ ap oc n 2 a 2 a 2 for small £, and 
WJ ap oc for large £. The binary regime C 1) 

is proportional to n 2 as expected. The n -1 / 2 dependence 
for £ 3> 1 clearly indicates the dominant contribution 
of phonon assisted transitions |TJ|,[Ï8|]. The increase in 
sound velocity s oc yfn with the condensate density n 
reduces the rate of phonon emission (and WJ ap ) [Q. 

As an example, we consider a sodium ion (Na^~) in a 
condensate of sodium àtoms at T = 0, with n ~ 10 14 
cm -3 and a — 52 a |2(J (a : bohr radius). For the 
uppermost bound state of the polarization potential, a v 
is well approximated by the atom-ion scattering length 
di ~ 2000 ao and assuming m/fj, « 1 0, we ob- 

tain f ~ 0.066, and neglecting n qo , WJ ap ~ 600 s _1 (see 
Fig. D). Under these conditions, roughly 600 àtoms will 
be captured by the ion, and another 600 will be emitted 
as phonons (excited out of the condensate) per second. 
In view of typical condensate lifetimes, this represents a 
sizable transfer mechanism that should be observable. 

By comparison, the capture rate into more deeply 
bound states will be much smaller. E.g., using the ap- 
proach of LeRoy and Bernstein |22| , we find for the sec- 
ond uppermost level of the pure polarization potential 
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where K 4 = 4^2nh 2 /n r(5/4)/[r(3/4)C* 4 1/4 ]. Because 
a v -\ < a v ~ di, the capture rate will be accompanied by 
atom-like excitations with Wy ap oc a 2 ,, and 
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For Na, C4 = 162.7 a.u. ||, and assuming m/^í w 1, 
we find that WJ ap is about ten times larger for the up- 
permost level than for the second uppermost level. 

The ratio is even larger for deeper levels: spontaneous 
capture is important mostly in the uppermost level. 

Let us now examine the escape rates out of the bound 
states v, Wy P and W^° wn . Since only the capture in the 
uppermost level is relevant, we consider only this level. 
At very low temperatures, k B T <C As, and there are no 
thermal phonons with sufficient energy to promote bound 
àtoms to the free condensate state, and non-equilibrium 
phonons emitted írom the capture process are moving far 
away from the location of the ion and can be neglected. 
Hence, because of energy conservation, transitions from 
bound levels to the continuum (into the condensate) are 
not allowed, and W^ p ~ 0. The collision of àtoms within 
the upper bound level (or condensate àtoms with the 
trapped ones) may result in the decay of the upper level 
to a lower level. The rate of this process (W^ own ) is in- 
versely proportional to the binding energy of the deeper 
level: W^ p /W^ OV!n ~ al/a 2 ^. Hence, M^ down is of the 
same order as W^|, and is at least one order of magni- 
tude smaller than W£ ap - 

As the number of trapped àtoms N v within the up- 
permost bound level v increases, the atom-atom repul- 
sive meanfield energy [p3[ will grow, and will eventually 
"push" the energy level up. Using a simple treatment 
based of the Gross-Pitaevskii equation, one can show that 
the binding energy of the uppermost level v, for N v 3> 1, 
is approximately given by e v (N v ) ~ (ma v /&jj,aN v ) 2 / 3 £ v . 
Although the number of àtoms trapped in the level v 
would apparently grow to extremely large numbers (since 
the decay rates out of it are much smaller), other effects, 
such as charge hopping or thermal fluctuations, will 
limit that maximum number. The most important one is 
related to thermal fluctuations. If their energy (~ k B T) 
is equal to |e,,(iV u )|, thermal equilibrium is reached, and 
as many àtoms get kicked out of the level as there are be- 
ing captured. In other words, W^ p will be equal to W^ ap . 
At that point, A^ nax can be found from e v (N™ ax ) ~ k B T, 
and Ç will have reached its maximal value £ ~ /ic/^T. 
For a BEC of Na àtoms with a v ~ R c ~ 2000ao and 
T ~ 100 nK, we get N™ ax ~ 600. In Fig. |, we illustrate 
the trajectory of the system as a function of £ for var- 
ious initial vàlues of a v , for both sodium and rubidium 
condensates. In all cases, as £ grows from its initial to 
its final equilibrium value, the capture rate W£ ap passes 



through a maximum before reaching its final value: then 
the rates "in" and "out" of the uppermost level are equal, 
and the system has reached thermal equilibrium. 

Làsers can also be used to control coherent stimulated 
transitions of free condensate àtoms into a specific bound 
state v of the molecular ion. Two off-resonant làser fields 
with a frequency difference matching the transition en- 
ergy into a specific bound state v cause a stimulated two- 
photon transition from the condensate into the bound 
state of molecular ion (see Fig. |ï|b). We assume that two 
làser beams are co-propagating, hence negligible momen- 
tum is transfered to the àtoms by the làsers. In the limit 
of sufficiently slow excitation, the atom-laser ficld intcr- 
action is described by H int = hCl J d 3 r^ v (r)*^ (r)blb + 
h.c, where Í7 is the Rabi-frequency of the two-photon 
transition proportional to the làser intensity, and bo,b v 
are bose annihilation operators for the condensate and 
bound àtoms, respectively. The resulting dynamics cor- 
responds approximately to that of a two-component con- 
densate with effcctivc coupling. The number of àtoms 
stimulated into the bound state can be estimated by solv- 
ing the equations of motion for bo,b v . For a sufficently 
short interaction time r, assuming an undepleted conden- 
sate (%bo) = N, and neglecting spontaneous relaxation, 
we obtain for N v = 



N v (t) ~ |Q|V/„(0) = 87ro;|íï|VJV/V • 



(10) 



The rate of stimulation W$" is therefore proportional to 
the làser power with a characteristic frequency scale given 
by ~ (íly/a^N/V), which can be easily made larger 
than spontaneous capture rate WJ ap . 

In practice, very fast excitation by làsers might result 
in heating of the condensate. The heating rate W^ eat is 
given by the transition probability from the bound state 
into continuum (Fig. |ï]c). This rate is also proportional 
to the làser power and can be estimated using Fermi's 
Golden rule. Taking as a characteristic transition 

time, we find that heating is negligible in dilute conden- 
sate (Na 3 /V < 1): W^ heat < W.f provided that 



W? < h/m{N/V) 2/z . 



(11) 



For realistic parameters, such as those of a sodium con- 
densate used above, the stimulated rate W** can be in 
the range of lOOKHz - 1 MHz without substantial heating 
of the condensate. For weakly bound states, such rates 
are achievable with modest làser powers used, e.g., in 
Bragg spectroscopy experiments pij . It is important to 
note that as the number of bound àtoms increases, atom- 
atom interaction can shift the bound state energy out of 
resonance, thereby inhibiting coherent transitions. This 
is directly analogous to the recently proposed "dipole blo- 
cade" in Rydberg àtoms |25| , and can be potentially used 
for non-trivial manipulation of quantum states in a BEC. 

In conclusion, we predicted the possiblity of forming 
large molecular ions in a BEC, analyzed the dynamics 
of their spontaneous growth and proposed a techniquc 
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for their fast and coherent generation in specific states. 
Such mesoscopic molecular ions could be observed, e.g., 
in experiments involving weak ionization of BECs. Of 
furthcr interest are experiments in which charged impu- 
rities are introduced via controlled processes involving 
either a combination of ion and atom traps, or using sur- 
face microtraps. The present work suggests that in such 
situations, iònic impurities can be used as a new tool 
for the accurate manipulation of condensates, including 
their quantum state. Finally we note that one could ma- 
nipulate, using external elèctric and magnètic fields, the 
position of the atòmic cloud trapped around a charged 
particle. 
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FIG. 1. Diagrams of atom capture by an ion. In (a), 
the spontaneous capture in level v is followed phonon emis- 
sion (with corresponding rates). In (b), the stimulated rate 
from the condensate c to the bound level v, also may produce 
hcating, i.e. population of noncondensate àtoms nc, in (c). 
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FIG. 2. WJ ap as a function of Ç for various vàlues of a v 
corresponding to the states 2 £jJ~, and 2 Ej of NaJ |J, as well 
as a v ~ R c for both Na and Rb. The equilibrium points for 
Na and Rb at 10 nK are shown, as well as the two limiting 
behaviours. 
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